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ABSTRACT 

The Magellanic Mopra Assessment (MAGMA) is a high angular resolution ^^CO{J = 
1 — >■ 0) mapping survey of giant molecular clouds (GMCs) in the Large and Small 
Magellanic Clouds using the Mopra Telescope. Here we report on the basic physical 
properties of 125 GMCs in the Large Magellanic Cloud (LMC) that have been sur- 
veyed to date. The observed clouds exhibit scaling relations that are similar to those 
determined for Galactic GMCs, although LMC clouds have narrower linewidths and 
lower CO luminosities than Galactic clouds of a similar size. The average mass sur- 
face density of the LMC clouds is 50 Mopc~^, approximately half that of GMCs in 
the inner Milky Way. We compare the properties of GMCs with and without signs 
of massive star formation, finding that non-star-forming GMCs have lower peak CO 
brightness than star- forming GMCs. The average CO-to-H2 conversion factor, Xqo, 
of non-star-forming GMCs is also ~ 50 per cent larger than for star-forming GMCs. 
We compare the properties of GMCs with estimates for local interstellar conditions: 
specifically, we investigate the H i column density, radiation field, stellar mass sur- 
face density and the external pressure. Very few cloud properties demonstrate a clear 
dependence on the environment; the exceptions are significant positive correlations 
between i) the Hi column density and the GMC velocity dispersion, ii) the stellar 
mass surface density and the average peak CO brightness, and iii) the stellar mass 
surface density and the CO surface brightness. The molecular mass surface density 
of GMCs without signs of massive star formation shows no dependence on the local 
radiation field, which is inconsi s tent w ith the photoionization-regulated star forma- 
tion theory proposed bv lMcKed ^M)- We find some evidence that the mass surface 



de nsity of the MAG MA clouds increases with the interstellar pressure, as proposed 
by Hi megreeni but the detailed predictions of this model are not fulfilled once 



estimates for the local radiation field, metallicity and GMC envelope mass are taken 
into account. 
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1 INTRODUCTION 



In th e Milky Way, molecular gas is mostly located in giant molecular clouds (GMCs) with masses M > 10^ Mq I Solomon et al 
19871 . henceforth S87). Understanding the physical properties of the GMCs is important because these clouds are the primary 
sites of star formation: the formation of GMCs and the transformation of molecular gas into stars are key processes in the 
life cycle of galaxies. Models of galactic evolution typically assume that GMCs are sufflciently similar across different galactic 
environments that a galaxy's star formatio n rate can be parameterised as the product of the CMC formation rate and the star 
formation efficiency of molecular gas (e.g. Ballesteros-Paredes fc HartmanrJ 2007 : Blitz fc Rosolowskv 2006h . This approach 
was initially justified by studies of Galactic molecular clouds, which found that the basic physical properties of GMCs in the 
Milky Way's disc obeyed well-defined scaling relations, often referred to as "Larson's laws" (e.g. S87, Larsorj|l98ll : Hever et al 



2001f ). More recently ^ considerable effort has been devoted to determ i ning whether GM Cs in other galaxies also follow the 
Larson relations (e.g. Rosolowskv et ai1l2003 : Rosolowskv fc Blitj|2005 : Rosolowskvll2007 ) . since empirical evidence that CMC 
properties are uniform - or at least exhibit well-behaved correlations with a parameter such as metallicity or pressure - would 
provide valuable information for developing models of star formation and galaxy evolution through cosmic time. 



As well as furnishing galaxy evolution models with empirical inputs, studies of extragalactic CMC populations aspire to 
resolve long-standing questions about the physical processes that are important for the formation and evolution of molecular 
clouds: are GMCs quasi-equilibrium structures, for example, or transient features in the turbulent interstellar medium? Do 
all GMCs form stars, and if not, why not? What is the physical origin of Larson's scaling relations? Although a nu mber of 
different theories to explain molecular cloud properties and the Larson relations have been proposed (e.g. M89, E89, IChieze 
19871 : iFleckl 1988 1. there are few extragalactic CMC samples that are comparable to the S87 catalogue, which contains 273 
clouds in the Galactic disc between longitudes 8° and 90°, and with radial velocities between -100 and 200 kms~^. The survey 
of ^^CO(J = 1 -> 0) emission in the LMC by NA NTEN (henceforth "the NANTEN survey") provided the first complete 
inventory of GMCs in any galaxy I Fukui et al. 2008h. but did not resolve molecular cloud structures smaller than ~40 pc (we 



adopt 50.1 kpc for the distance to the LMC, e.g. AlvesI 2004 ) . Thorough testing of the different molecular cloud models will 



require deep, unbiased wide-field surveys of molecular clouds at high angular resolution across a range of interstellar condi- 
tions. Extensive surveys of this kind are only just feasible with current instrumentation, and hence the number of molecular 
cloud samples that can be used to falsify molecular cloud models remains frustratingly small. 



To date, studies of the CO emission in nearby galaxies have concluded that extragalactic G MCs are alike. For a sample of 
~ 70 resolved GMCs located in five galaxies (M31, M33, ICIO, and the Magellanic Clouds). [Blitz et al.l l|2007l ') found that 
extragalactic GMCs not only follow the Galactic Larson relation s , but also that different galaxies have similar GMC mass 
distributions. Similar conclusions were reached by Bolatto et al. ( 20081 . henceforth B08) using a sample of ~ 100 resolved 
GMCs in twelve galaxies, although these authors noted that molecular clouds in the SMC tend to have low CO luminosities 
and narrow linewidths compared to GMCs of a similar size in other galaxies. By comparing tracers of star formation and 
neutral gas on ~ 1 kpc sca le s for galaxies in Th e If i Nearby Galaxy Survey (THINGS, which does not include the Magellanic 



Clouds, IWaher et al.ll2008l l. iLerov et al 



( 200^ found that the star-forming efficiency of molecular gas (defined as the star 
formation rate surface density per unit molecular gas surface density SFEn2 = Sspr/Eh2) is approximately constant in 
normal spiral galaxies, SFEn^ = 5.25 ± 2.5 x lO"'^" yr~^. As noted by the authors, this result could arise if the star-forming 
efficiency of an individual GMC is d etermined by its intrinsic properties, and if the properties of GMCs are independent 
of their interstellar environment (e.g. Krumholz fc McKeel 20051 ). While the existing observational evidence has so far been 
interpreted in favour of uniform GMC properties, a dependence of GMC properties on the local interstellar environment is 
by no means ruled out. A constant SFE^^ on kiloparsec scales indicates that the properties of GMC ensembl es are alike 



on tho se scales; whether this conclusion can be applied to individual GMCs is far less certain. Neither B08 nor iBlitz et al 



( 20071 ) pursued the origin of the scatter in the extragalactic Larson relations that they observed, moreover, even though the 
mean GMC mass surface density for the galaxies in their respective samples varies by more than an o rder of magnitude , and 



the mass surface densities of the individual GMCs varies between ~ 10 and 1000 M0 pc ^ (see also Hever et al. 20081 . for 



evidence that the mass surface density of Milky Way clouds is not constant). A resolved survey of a large number (> 100) of 
GMCs located in a single nearby galaxy therefore remains valuable, since it eliminates the uncertainties inherent in combining 
heterogeneous datasets and provides a sample that is large enough to investigate both the average properties and scaling 
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relations of an extragalactic GMC population, as well as the dispersion around overall trends and average quantities. 



In this paper, we report on some initial results from the Magellanic Mopra Assessment (MAGMA), an ongoing, high-resolution 
survey of the 12C0(J =1^0) emission from molecular clouds in the Magellanic Clouds using the Mopra Telescope. Here 
we present results from the Large Magellanic Cloud (LMC) only; a description of the molecular clouds surveyed by MAGMA 
in the Small Magellanic Cloud (SMC) has been presented elsewhere (MuUer et al., accepted). While the ^^CO(J = 1 — >■ 0) 
emission from molecular gas in the LMC has been the target of extensive mapping with the NANTEN telescope and Swedish- 
ESO Submill imetre Telescope (SEST), neither pro ject obtained observations that were ideal for studying the Larson relations 
in the LMC I Fukui et al. 2008 : Israel et al. 2003h. The spatial resolution of the NANTEN survey is comparable to the size 



of a typical Milky Way GMC (~ 50 pc, e.g. Blitz 1993l l: ideally, we would like to resolve structures on smaller spatial scales 



in order to include less massive GMCs in our analysis. Resolved observations are crucial, moreover, for accurate estimates 
of derived GMC quantities such as virial mass and mass surface density. The SEST Key Programme CO in the Magellanic 
Clouds mapped molecular clouds in the LMC with comparable spatial resolution as MAGMA (~ 10 pc), but was strongly 
biased towards regions associated with well-known sites of active star formation. An analysis of the striking molecular cloud 
complex s i tuate d so uth of the 30 Dorad us star-forming complex using the MAGMA data has already been presented by 
Ott et al.l (|2008l l and lPineda et al.l \20od ): in this paper, we turn our attention to the general LMC cloud population. 



This paper is structured as follows: in Section [21 we summarise the MAGMA observing strategy and our data reduction 
procedure, and also describe the ancillary data that we have used in our analysis. Section [3] outlines the approach that we 
have used to identify GMCs and to measure their physical properties. In Section[4l we compare the properties of GMCs with 
and without star formation. Scaling relations between the cloud properties are discussed in Section[5l while Section[6]presents 
a comparison between the intrinsic physical properties of the GMCs and properties of the local interstellar environment. In 
SectionjTl we discuss whether our results are consistent with i) the photoionization-regulated theory of star formation proposed 
by McKeel ( 1989I . hencefort h M89) and ii) a d ominant role for interstellar gas pressure in the determination of molecular cloud 
properties, as suggested bv lElmegreen j 19891 . henceforth E89). We conclude with a summary of our key results in Section [S] 



2 DATA 



2.1 The MAGMA LMC survey: observations and data reduction 

MAGMA observations of the ^^CO(J =1-^0) emission from molecular clouds in the LMC are conducted at the Mopra Tele- 
scope, which is situated near Coonabarabran, Australia^ At 115 GHz, the Mopra Telescope has a full width at half-maximum 
(FWHM) beam size of 33 arcsec, corresponding to a spatial resolution of 8 pc at our assumed distance to the L MC. Due to the 
large angular size of the LMC's gas disc , and the small cov ering fraction of the ^^CO(J =1^0) emission (e.g. lMizuno et al.l 
200ll ). we use the NANTEN survey of IPukui et"al] (|2008t ) to select regions of bright ^^CO(J = 1 ^ 0) emission for high- 
resolution mapping. MAGMA observations target all molecular clouds with CO lumino sities Leo > 7000 Kkms~^ pc'^ and 



peak integrated intensities > 1 Kkms~^ in the NANTEN cloud catalogue of I Fukui et al.l (2008). In combination, these clouds 
contribute ~ 70 per cent of the LMC's total CO luminosity. An integrated intensity map of the CO emission in the LMC by 
MAGMA is shown in Fig. [1] with outlines that represent the survey's coverage. 



MAGMA observations are conducted in "on-the-fly" (OTF) raster-mapping mode. A grid of 5 x 5 arcmin^ OTF fields is 
placed over the region surrounding a molecular cloud target. The field centres are separated by 4.75 arcmin in RA and Dec; 
this 15 arcsec overlap ensures full coverage of the molecular cloud and facilitates mosaicking. In OTF mode, the telescope 
takes data continuously while scanning across the sky. Along each row of an OTF field, individual spectra are recorded ev- 
ery 14 arcsec, so that the telescope beam is oversampled in the scanning direction. The spacing between rows is 10 arcsec, 
also oversampling the beam. Each row is preceded by an off-source (emission-free) integration situated approximately 10 
to 20 arcmin from the field centre; an absolute location for the OFF spectra was selected for each molecular cloud target 
in order to ensure a consistent sky subtraction. To minimize scanning artefacts, each field in the final survey dataset is 
mapped twice: an initial pass is made with the telescope scanning in the RA direction, and a second pass then conducted 
in the orthogonal direction. The data presented in this paper were obtained during the southern hemisphere winters of 2005 
to 2008. Observations are complete for ~ 40 per cent of the clouds presented here; the remaining clouds have been scanned 
in the RA direction only. The full MAGMA LMC dataset will be presented in a subsequent paper (Wong et al, in preparation). 



^ The Mopra Telescope is managed by the Australia Telescope, which is funded by the Commonwealth of Australia for operation as a 
National Facility by the CSIRO. 
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Figure 1. A map of CO integrated intensity in the LMC by MAGMA (greyscale). The thin black Uncs indicate the coverage of the 
MAGMA observations included in this paper. The black cross indicates the position of 30 Doradus; in the text, we refer to the large 
survey region that encloses and extends south from 30 Dor as the "molecular ridge" . The black ellipse indicates where the stellar surface 
density, S*, is greater than 100 Mq pc~^; we refer to this region in the text as the "stellar bar". 



A system temperature measurement is obtained with an ambient load at the start of each OTP map and every 20 to 30 
minutes thereafter. In the interval between these measurements, the system temperature is monitored using a noise diode. 
Typical system temperatures for the survey observations are between 500 and 600 K; observing is abandoned when system 
temperatures exceed ~1000 K. Our intention is to obtain data with uniform sensitivity, so OTF fields with average system 
temperatures above ~ 850 K are re-observed. Between each OTF map, the pointing of the antenna is verified by observing 
the nearby SiO maser RDor. The pointing solution is updated (and re-verified) if the pointing error in azimuth or elevation 
is greater than 8 arcsec. Prior to correction, the pointing errors are typically below 10 arcsec. 



In 2005, the Mopra Telescope was equipped with a dual polari sation SIS receiver that produced 600 MHz instantaneous 
bandwidth for observing frequencies between 86 and 115 GHz jMoorev et al.l 119971 ). Th e correlator at that time could be 
configured for bandwidths between 4 and 256 MHz across the receiv er 's 600 M Hz band dW ilson ct al.| ll992l ). Our 2005 ob- 
servations targeted the molecular ridge region discussed by Ott et al. ( 20081 ) and Pineda et al. (2009) (see Fig. [TJ, with some 
data obtained for additional fields near RA 05hl6m, Dec -68dl0m (J2000) and RA 05h24m, Dec -69d40m (J2000) (J2000). 
For all 2005 observations, the correlator was configured with 1024 channels over a 64 MHz bandwidth centred on 115.16 MHz, 
which provided a velocity resolution of 0.16 kms"'^ per channel across a reliable velocity bandwidth of ~ 120 kms~^. As this 
is not quite sufficient to cover the total radial velocity range of the LMC's CO emission, the centre of the observing band was 
placed near the peak of CO spectrum obtained by NANTEN for the region being observed. In subsequent years, the data were 
recorded using the newly installed MMIC receiver and the University of New South Wales Digital Filter Bank (MOPS)El In 
the narrowband configuration used for the MAGMA survey observations, MOPS can simultaneously record dual polarisation 



^ The University of New South Wales Digital Filter Bank used for the observations with the Mopra Telescope was provided with support 
from the Australian Research Council. 
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data for up to sixteen 138 MHz windows situated within an 8 GHz band. Each 138 MHz window is divided into 4096 channels; 
for our survey observations at 115 GHz, this configuration provides a velocity resolution of 0.09 kms^^ per channel across 
the velocity range [90,410] kms~^. 



The Mopra beam has been described bv lLadd et al. I I2OO5I ). These authors identify three components - the 'main beam', and 
the inner and outer 'error beams' - that contribute to the antenna response. The presence of the se non-negligible er ror beams 
implies that the telescope efflciency, rj, will depend on source size. In the 2004 observing season, iLadd et al.l (|2005l ) derived a 
main-beam efficiency factor of rjmh = 0.42 for observations at 115 GHz, and an 'extended' beam efficiency of 77xb = 0.55 for 
sources that are comparable to the size of the inner error beam (~ 80 arcsec). For clouds in the MAGMA survey, we consider 
the extended beam efficiency factor to be more appropriate. 



To monitor the reliability of the flux calibration, we observed the standard source Orion KL (RA 05h35ml4.5s, Dec - 
05d22m29.56s (J2000)) once per observing session throughout our survey. In 2005, we measured an average peak brightness 
temperature of Orion KL of 55 ± 3 K (in TX units). According to the SEST documentationlfl the peak antenna temperature 
of Orion KL is ~ 7 1 K, corresponding to a main beam temperature of Tmb ~ 102 K for an efficiency factor of r^mb = 0.7 



( Johansson et al. 



199^ ). This suggests that data from the two telescopes can be placed on the same brightness temperature 
sc ale using a factor of rjxb = 0.54 ± 0.03 for the Mopra data, in excellent agreement with the telescope efficiency derived 



bv iLadd et al.l (|2005l ). In subsequent years, the average peak brightness temperature that we have measured for Orion KL 
has varied. In 2006, we measured an average peak brightness temperature of 35 ± 4 K, while our average measurements in 
2007 and 2008 were 43 ± 3 K and 50 ± 3 K respectively. For these data, we therefore used conversion factors of 0.35 (2006), 
0.43 (2007) and 0.49 (2008). To validate the fi nal flux scal e of ou r Mopra data, we c o mpiled published SE S T measurements 
of LM C m olecular clouds from the studies by Israel et al. ( 1993l l. Chin et al. (1997), Kutner et al. 1 1997 ). Johansson et al 



1 19981 ) and Hsr^Lelal] (|2003l ). We were able to compare Mopra and SEST measurements of the peak brightness and peak 
integrated intensity for 40 clouds. For both the peak brightness and the peak integrated intensity of these clouds, the average 
ratio between the SEST and Mopra measurements was 1.1, with a dispersion of ~ 20 per cent. 



The processing of the survey data involved four main steps. An initial correction to align the position and time stamp 
informat ion in the raw d ata flies was applied to data obtained in 2005, using the mopfix task of the MIRIAD software 
package (jSault et al.lll995l ). Bandpass calibration and baseline-fitting were performed using the AIPS++ livedata package. 
LiVEDATA determines the bandpass calibration for each row of an OTF map using the preceding OFF scan and then fits a 
user-specified polynomial to the spectral baseline. We chose to fit the spectral baselines with a first-order polynomial: data 
subcubes containing baselines showing higher order ripples were rejected. The spectra were combined to form a spectral line 
cube using the AIPSH — h GRIDZILLA package. Each [x,y,v) cell within the final cube is sampled multiple times by the OTF 
scans: GRIDZILLA averages the spectra contributing to the emission within each cell according to a convolution kernel, beam 
profile and weighting scheme specified by the user. We chose to grid the data with a cell size of 9 x 9 arcsec^, aggregating the 
data from both polarisations and both scanning directions. The data were weighted by the inverse of the system temperature 
measurements. We used a truncated Gaussian convolution kernel with a FWHM of 1 arcmin and a cutoff radius of 30 arcsec, 
producing an output data subcube with effective angular resolution of 45 arcsec. Although GRIDZILLA can process data from 
contiguous OTF fields, it cannot process data recorded with different correlator configurations simultaneously. The data sub- 
cubes from different years were thus processed separately, and were converted to Tmb units using the annual extended beam 
efficiency factors derived from our Orion KL observations. The MIRIAD task imcomb, which weights the input data by the 
inverse of the RMS noise, was then used to combine data from regions that were observed over multiple observing seasons. 
Finally, in order to increase the signal-to- noise (S/N), we binned the MAGMA data subcubes to a velocity resolution of 
0.53 kms~^ and smoothed them to an angular resolution of 1 arcmin. The average RMS noise in these subcubes is 0.24 K per 
0.53 kms^^channel. An example of the integrated ^^CO(J = 1 — >■ 0) emission for one such MAGMA subcube is shown in Fig.[2l 



2.2 Ancillary data 

2.2.1 Hi data 



To trace the atomic gas in the LMC, we use the H i rnap miblished bv iKim et all |2003t ). which combines data from the 



Australia Telescope Compact Array (ATCA Kim et al. 19981 ) and the Parkes single-dish telescope ( Stavelev-Smith et al 



2OO3I ). The angular resolution of the H i data is 1 arcmin, well- matched to the reduced MAGMA CO data. The H i data cube 



has a velocity resolution of 1.65 kms ^, with a column density sensitivity of 7.2 x 10 cm per channel. We construct an 



' |http: //www. apex-telescope ■ org/sest/html/telescope-calibration/calib-sources/orlonkl .html| 
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Figure 2. A ^^CO(J = 1—^0) integrated intensity map of an example GMC in the MAGMA cloud list. The contour spacing is 2 kms~^, 
with the lowest contour at Iqq = 1.5 kms~^. The ellipse represents the fit derived by CPROPS, from which the major and minor axis 
and position angle of the GMC were taken. As explained in Section|3l we are interested to identify GMCs in the MAGMA data subcubes, 
rather than structure on smaller scales corresponding to the spatial resolution of the survey. The black circle in the lower right corner 
represents the Mopra beam. 



Hi integrated intensity emission map of the LMC, /hi, by integrating the Hi data cube over the heliocentric velocity range 
196 to 353 kms"\ We assume that the Hi emission is optically thin everywhere, and derive a map of the LMC's Hi column 
density, A''(HI), according to 



Af(HI)[cm-2] = 1.82 x 10^*/Hi[Kkms"^]. 



(1) 



Our estimate for A'^(Hi) is likely to be a lovyer hmit sinc e the Hi emission i n the LMC may have sign i ficant optical depth , 
especially along the sightlines towards molecular clouds ( Dickey et al. 1994 : Marx-Zimmer et al. 2000l : Bernard et al. 
see also Section [7.2.2p . 



2.2.2 Stellar mass surface density 



T o trace the mass di stribution within the LMC's stellar disc, we use the stellar mass surface density map presented in fig. Ic 
of I Yang et al.l (120071 ). The map is based on number counts of red giant branch (RGB) and asy mptotic giant branch ( AGB) 
stars, selected by their colours from the Two Micron AH Sky Survey Point Source Catalogue ( Skrutskie et al. 20061 ) . The 
star counts are binned into 40 pc x 40 pc pixels, and then convolved with a Gaussian smoothing kernel with a = 100 pc. 
The resulting map is nor malized to a me asure of absolute stellar mass surface density by adopting a total stellar mass for 
the LMC of 2 X 10^ M© ( Kim et al. 19981 ). The resolution of the stellar surface density map is considerably coarser than the 



angular resolution of our MAGMA data subcubes, so it is possible that the average stellar mass surface density of smaller 
GMCs is underestimated due to beam dilution. RGB and AGB stars are relatively old populations, however, so their spatial 
distribution is likely to be smooth. In p articular, we do not expect them to be strongly clustered in the vicinity of molecular 



clouds (e.g. iNikolaev fc Weinberg||200i 



2.2.3 Interstellar radiation field 

To estimate the interstellar r adiation field at the lo cations of molecular clouds within the LMC, we use the dust temperature 
(Tfj) map prese nted in fig. 7 of BerMJx]_et_al. ( 200^ . The map is derived u sing the ratio of the IRIS 100 pm and Spitzer 160 pim 
emission maps ( Miville-Deschenes fc Lagache 2005 : Meixner et al. 20061 ) ■ assuming that dust emission can be modelled as a 
grey body: 



(2) 
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with P — 2 at far-infrared wavelengthsR lBernard et alj (|2008l ) find tliat dust temperatures in the LMC vary from 12 K up to 
34.7 K, with an average value of 18.3 K. This is significantly colder than previous determinatio ns, especially f o r estim ates that 
attempted to constrain the dust temperature using the IRAS 60/im flux density. As noted by Bernard et al 1 2008h . emission 
at 60/im is highly contaminated by out-of-equilibrium emission from very small grains (VSGs), and this is especially true in 
the LMC, due to the presence of excess 70^m emission. Temperatures derived from the IRAS 60/100/im flux density ratio 
may therefore be strongly overestimated. 



For dust grains in t hermal equilibrium, the strength of the radiation field. Go, is related to the dust temperature by Go oc 



4+/3 



(e.g. Leaueujboosi ). The average strength of Go across the entire LMC i s thus only a factor of ~ 1.3 greater than in the solar 
neighbourhood. Go,©, where the dust temperature is 17.5 K (with (3 — 2, Boulanger et al. 1996h . At the locations observed by 
MAGMA, Go/Go,0 varies between 0.5 and 58.8, with a median value of 1.7. For Go/Go,0 averaged over the projected cloud 
areas, we obtain values in the range 0.5 to 8.0. Our derived value for the strength of the radiation field near the well-known 
LMC star- formin g region N159W (Go/Go n = 8.4) is very low compared to the value obtained by previous estimates (e.g. 



Israel et al 



1996; 



Pineda et allboosl . I2OO9I ) . Part of the discrepancy may be due to beam dilution since the resolution of the 



Bernard et al.. (12008 ) dust temperature map is 4 arcmin, which corresponds to a spatial scale of 60 pc in the LMC. However, 
we also note that our method for determining Go/Gq,© is extremely sensitive to the as sumed dust temperat ure, which is colder 
than reported by previous analyses. At the locations of the molecular clouds in the [Bernard et al. I I2OO8I ) dust temperature 
map, the mean (maximum) formal error on Td is 2 per cent (12 per cent). These errors do not include potential variations 
in the value of j3, however, so the true uncertainty is likely to be greater than this. Despite the uncertainties regarding the 
absolute value of Go/Go,q, we emphasise that our analysis concerns the relationship between CMC properties and the relative 
strength of the radiation field, and our conclus ions do not rely on an absolute calibration of Go/Gq,©. Finally, we note that 
while the resolution of the Bernard et al. ( 20081 ) dust temperature map is coarser than the angular resolution of our MAGMA 
CO data, it is still likely to provide a reasonable estimate for the average dust temperature within individual GMCs with 
R ^ 30 pc. For smaller clouds, the average dust temperature and average radiation field may be underestimated due to beam 
dilution. 



3 MEASURING GMC PROPERTIES 



To identify GMCs in the M AGMA data subcubes and measure their properties, we have used the algorithms presented in 
Rosolowskv &L Lerov l|2006l . implemented in IDL as part of the CPROPS package). CPROPS uses a dilated mask technique 



to isolate regions of significant emission within spectral line cubes, and a modified watershed algorithm to assign the emission 
into individual clouds. Moments of the emission along the spatial and spectral axes are used to determine the size, linewidth 
and flux of the clouds, and optional corrections for the finite sensitivity and instrumen tal resolution may be appli ed to the 
measured cloud properties. Each step of the CPROPS method is described in detail by Rosolowskv fc Lerov 1 2006[ ). 



Regions of significant emission within the MAGMA data subcubes are initially identified by finding pixels with emission greater 
than a threshold of 4cri{MS across two contiguous velocity channels. The mask around these core regions is then expanded to in- 
clude all the pixels connected to the core with emission greater l.baRMS across at least two consecutive channels. We explored 
a range of values for the threshold and edge parameters in the masking process, and found that these values distinguished 
credible emission regions (i.e. the mask did not expand excessively into the noise) and also yielded reliable measurements 
for the properties of faint clouds. The emission identified with an isolated cloud in the MAGMA dataset is illustrated in Fig. [S] 



Once regions of significant emission have been identified, CPROPS assigns the emission to individual cloud structures. To 
generate the preliminary list of GMC can didates, we used the default parameters for the identification of GMCs that are rec- 
ommended in iRosolowskv fc Lerovl (|2006l l. In this case, the parameters of the decomposition are motivated by the observed 
physical properties of Galactic GMCs: spatial sizes greater than ~ 10 pc, linewidths of several kms~^, and brightness tem- 
peratures less than ~ 10 K. We adopt this approach because our goal is to describe the properties of G MCs in the LMC and 



to inv estigate how these properties might differ to the properties of GMCs in other galaxies. As noted bv IRosolowskv fc Lerov 
( 2006l l. GMCs contain structure across a wide range of size scales, so identifying the clumpy substructure - with a size scale of 
~ 1 pc and typical linewidth of ~ 1 km s~^ - within the clouds in our spectral line cubes would require different decomposition 
parameters than the ones that we have used. While the properties and scaling relations of this substructure is an important 
topic for investigation, we defer this analysis to a future paper. A previous anal ysis of the MAGMA data in the molecular ridge 



region employed a difi'erent decomposition algorithm (i.e. GAUSSCLUMPS, 



Pineda et al 



2OO9I ), 



which identified structures 



■* 'IRIS' is an abbreviation for Improved Reprocessing of the IRAS Survey. As described in iMiville-Deschenes fc Lagachj l|2005l '). IRIS 
images have improved zodiacal light subtraction, absolute calibration and scanning stripe suppression than the IRAS Sky Survey images. 
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Figure 3. Channel maps of the CO emission from the CMC in Fig. [2] (greyscale). The black contour indicates the emission region that 
CPROPS identifies as belonging to the cloud. The velocity axis of the MAGMA data subcube and CPROPS assignment cube has been 
binned to a channel width of 2.1 kms^^ for illustration only. The black circle in the lower right corner of each panel represents the 
Mopra beam. 



that are typically smaller (with radii between 6 and 20 pc) than the GMCs described here. 



Our initial list of GMC candidates contained 237 objects. After cross-checking the CPROPS cloud identifications against the 
NANTEN ^^CO(J = 1 — >■ 0) data cube and the MAGMA data subcubes to verify that the identified objects were genuine, we 
rejected 54 features that were clearly noise peaks or map edge artefacts. To ensure that the properties of clouds in our final 
list are reliable, we impose a signal-to-noise threshold S/N ^ 5 and reject objects with measurement errors for the radius 
and velocity dispersion that are greater than 20 per cent. We further require that CPROPS is able to apply the corrections 
for sensitivity and instrumental resolution successfully (i.e. these corrections do not result in undefined values). The resulting 
sample of GMCs, which we refer to as the "MAGMA cloud list" , contains 125 clouds. To verify the results of our analysis, 
we also define a "high quality" GMC sample, which contains 57 clouds with S/N J5 9 and for which the uncertainties in the 
radius and velocity dispersion measurements are less than 15 per cent. 



Cloud properties follow the standard CPROPS definitions. The cloud radius is defined as i? = 1.91cr_R, where an is the 
geometric mean of the second moments of the emission along the cloud's major and minor axes. The velocity dispersion, a^, 
is the second moment of the emission distributi on along the velocity axis, which for a Gaussian line profile is related to the 
FWHM linewidth, Av, by Av = \/8 In 2 CTv (e.g. Rosolowskv &: Lerov 20061 ) . The CO luminosity of the cloud, Leo, is simply 
the emission inside the cloud integrated over position and velocity, i.e. 



Lco(K kms 



-1 2, 
pc ) 



D- 



( - V 

V180 X 3600/ 



(3) 



,180 X 3600. 

where D is the distance to the LMC in parsecs, 5x and Sy are the spatial dimensions of a pixel in arcseconds, and 5v is the 
width of one channel in kms~^. The virial mass is give n by Mvir = 1040O-V-R M©, which assumes that molecular clouds are 
spherical with p oc density profiles ( MacLaren et al.1ll988h . CPROPS estimat es the error associated with a cloud property 
measurement using a bootstrapping method, which is described in section 2.5 of Rosolowskv fc Lerov I 2006h . 



In addition to the basic properties reported by CPROPS, we define the CO surface brightess as the total CO luminosity of 
a cloud divided by its area, Jco[Kkms~'^] = Lco/i^R^ ■ The molecular mass surface density, Eh2, is defined as the virial 
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mass divided by the cloud area, Ehj [Mq pc-2] = M^i./nR^. The CO-to-Ha conversion factor, Xco, may be expressed as the 
ratio between these quantities: Xco [cm~^ (Kkms~^) ^] = [Mq pc~'^]/(2.2/co [ Kkms~^]). We note that our definitions 
of Xco and assume that GMCs manage to achieve dynamic equilibrium: if the degree of virialisation is not constant for 
LMC molecular clouds, then variations in Xco and Ehj may instead reflect differences between the dynamical state of the 
clouds. We estimate the uncertainties in Xco, -'^co and using standard error propagation rules. 



As emphasized bv lRosolowskv &: Lerov l|2006l l. the resolution and sensitivity of a dataset influence the derived cloud proper- 



ties. In order to reduce these observational biases, they recommend extrapolating the cloud property measurements to values 
that would be expected in the limiting case of perfect sensitivity (i.e. a brightness temperature threshold of K), and cor- 
recting for finite resolution in the spatial and spectral domains by deconvolving the telescope beam from the measured cloud 
size, and deconvolving the channel width from the me asured linewidth. The procedures that CPROPS uses to apply these 
corrections are described in Rosolowskv Sz Lerov ( 2006I . see especially fig. 2 and sections 2.2 and 2.3). For the analysis in this 



paper, we use cloud property measurements that have been corrected for resolution and sensitivity bias. 

The GMCs observed by MAGMA in the LMC have radii ranging between 13 and 160 pc, velocity dispersions between 1.0 and 
6.1 kms~^, peak CO brightnesses between 1.2 and 7.1 K, CO luminosities between 10^'^ and 10^'^ K kms^^ pc^, and virial 
masses between 10*'^ and 10^'* Mq. The clouds tend to be elongated, with a median axis ratio of 1.7. A detailed analysis of 
the cloud property distributions and their uncertainties will be presented once the MAGMA LMC survey is complete (Wong 
et al, in preparation). 



4 PHYSICAL PROPERTIES OF NON-STAR-FORMING GMCS 



Molecular clouds have traditionally been modelled as quasi-equilibrium structures, but recent theoretical work has also begun 
to explore whether molecular clouds might form and disp erse more rapidly as a consequence of la rge-scale dynamical events 
in the ISM, such as turbulent fiows or cloud collisions (e.g. Bergin et al. 2004 : Tasker fc Tan 2008*) . Observational constraints 
on the physical properties of recently formed GMCs would be a useful contribution to the debate about GMC lifetimes, but it 
remains unclear whether observations of ^^CO(J = 1 — >■ 0) emission alone can distinguish between younger and more evolved 
GMCs. While it is plausible that some potential characteristics of newly formed GMCs, such as colder gas temperatures, 
stronger bulk motions or sparser CO filling factors, would have observational signatures, it is also possible that the onset of 
widespread ^^CO(J = 1 — >■ 0) emission might occur late in the cloud formation process, or that physical conditions in the 
CO-emitting regions of GMCs are relatively uniform and therefore insensitive to a cloud's evolutionary state. In this section, 
we investigate whether there are significant differences between the properties of young GMC candidates and other GMCs in 
the MAGMA LMC cloud list. 



We co nstructed a sample of young GMC candidates using th e evolutionary clas si fication scheme designed bv lKawamura et al 
( 2009^ for the 272 GMCs in the NANTEN LMC catalogue (jPukui et aklboOSl ). iKawamura et all (|2009h classified GMCs on 
the basis of their association with Hii regions and young stellar clusters, finding 72 Type I GMCs that show no association 
with massive star-forming phenomena. An important assumption behind this approach is that all GMCs eventually form 
stars and are finally dissipated by their stellar offspring; in this scenario, GMCs without signs of massive star formation may 
be considered young. To date, MAGMA has observed 30 Type I NANTEN GMCs, but there are only 17 MAGMA clouds 
associated with these Type I GMCs that satisfy our criteria for inclusion in the MAGMA cloud list: henceforth, we refer 
to these 17 MAGMA clouds as the "young GMC sample" or the "non-star-forming GMCs". The 17 clouds in the young 
GMC sample correspond to 16 NANTEN GMCs: 15 of the 17 MAGMA clouds demonstrate a one-to-one correspondence 
with a Type I GMC in the NANTEN catalogue, while the other NANTEN GMC divides into two clouds at MAGMA's finer 
resolution. CPROPS identifies 10 more objects in the MAGMA data subcubes that are associated with a further 8 NANTEN 
Type I GMCs, but these objects either have S/N < 5 and/or measurement errors that exceed 20 per cent, and are therefore 
excluded from the MAGMA cloud list. CPROPS does not identify any significant emission in the MAGMA data subcubes 
for the remaining six NANTEN Type I GMCs that we have observed to date, although we note that these are all regions of 
the MAGMA survey tha t have only been scann ed once. We do not attempt to re-classify the MAGMA clouds according to 
the criteria developed by Kawamura et al. ( 20091 ) . but simply ascribe the evolutionary classification of the NANTEN GMC to 
any MAGMA cloud that is coincident in space and velocity. This classification should be reliable for clouds in the young GMC 
sample since their projected areas are always smaller than, and contained within, the corresponding NANTEN GMC boundary. 



To verify that differences between the properties of star-forming and non-star-forming GMCs are not due solely to variations 
in the cloud size, we constructed a control sample with the same size distribution as the young GMC sample by matching 
each of the 17 young GMC candidates with three star-forming clouds in the MAGMA cloud list of a similar radius (i.e. 
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Table 1. Average physical properties of the MAGMA elouds and results of the KS tests for the young GMC sample. Columns 2 to 
4 list the median and median absolute deviation of th properties of clouds in the MAGMA cloud list, the young GMC sample and 
the control sample respectively. Column 5 lists the median P value obtained in the error trials, and column 6 lists the standard 
deviation of the P values in the trials (see text). 



Cloud Property 


MAGMA 


Young GMCs 


Control 




a(P) 




28±8 


30it8 


30it7 


0.96 


0.05 


/T", - ( LrTYl G ^\ 
V \ JVlii o 1 


2.3it0.5 


2.4±0.7 


2.1±0.5 


0.54 


0.18 




l.litO.l 


l.OiO.l 


l.litO.l 


0.007 


0.007 




2.2it0.4 


1.7it0.2 


2.3±0.4 


0.003 


0.003 


Leo (10'' K kms~i pc^) 


1.2±0.6 


1.2±0.5 


1.2±0.5 


0.89 


0.11 


Mvir (10^ Mo) 


1.5±0.9 


1.7±1.0 


1.5±0.8 


0.78 


0.18 


Axis Ratio, F 


1.7±0.4 


1.5±0.2 


1.8±0.4 


0.25 


0.10 


Sh^ (Mq pc-2) 


55±22 


73±21 


52±23 


0.13 


0.10 


Ico (K kms-1) 


4.8±1.5 


4.3±0.8 


4.6±1.3 


0.78 


0.15 


Xco (cm-2 (Kkms-i)"^) 


4.7±1.6 


6.9±2.3 


5.4±2.0 


0.06 


0.05 


Vlsr (kms~l) 


255±22 


248±18 


256±23 


0.18 


0.00 


-Rgal (kpc) 


l.SitO.g 


2.3±0.9 


1.8±0.8 


0.43 


0.00 


Af(Hl) (xlO^i cm-2) 


2.7±0.7 


2.7±0.4 


2.7±0.7 


0.66 


0.18 


S. (M0PC-2) 


48±24 


50±16 


42±19 


0.25 


0.10 


Go/Go, 


1.7±0.5 


1.3±0.5 


2.0±0.6 


0.003 


0.003 


(Tpk(HI)> (K) 


68±10 


51±10 


68±8 


0.0001 


0.0001 


Ph/fcB (10* K cm-3) 


6.4±2.5 


6.4±0.9 


5.54±2.1 


0.54 


0.18 



the control sample contains 51 clouds). iKawamura et al. found that young GMCs in the NANTEN catalogue tend to 



be smaller than star-forming GMCs, but here we wish to determine whether there are differences between star-forming and 
non-star-forming GMCs that persist even after variations in cloud size are suppressed. The average discrepancy between the 
radius of a young GMC candidates and their three corresponding control clouds was 2 per cent, with a maximum discrepancy 
of 8 per cent. The average properties of GMCs in the MAGMA cloud list, the young GMC sample and the control sample are 
listed in Table [1] 

To test whether young GMCs have distinct physical properties, we conducted Kolmogorov-Smirnov (KS) tests between the 
young GMC sample and the control sample. The KS test is a non-parametric test that compares the cumulative distribution 
functions of two samples in order to determine whether there is a statistically significant difference between the two popu- 
lations. The test is reliable when the effective number of data points, A^e = n'i^N2 ' greater than 4, where A'^i and N2 are 
the number of data points in the first and second samples. The result of the KS test can be expressed as a probability, P, 
that the sample distributions are drawn from the same underlying distribution. As the traditional KS test does not account 
for uncertainties in the measurements of the cloud properties, we performed 1000 trials of each KS test. For each trial, the 
value of each cloud property measurement was displaced by kAx, where Ax is the absolute uncertainty in the cloud property 
measurement and A; is a uniformly distributed random number between -1 and 1. A summary of the KS test results is shown 
in Table [1] If the results of the error trials are narrowly distributed around zero - i.e. if (P) ^ 0.05 and (t(P) ^ 0.05, where 
cr represents the standard deviation - we consider that there is statistically significant evidence against the null hypothesis 
that the cloud samples are drawn from the same underlying population. 

Table[l]shows that the properties of both the young GMC and the control sample show a large dispersion around their average 
value, and that differences between the samples are not always obvious from measures of central tendency. The results of the 
KS tests indicate that the null hypothesis cannot be rejected for many of the intrinsic cloud properties. Both the maximum 
peak CO brightness within a GMC (Tmax) and the average peak CO brightness ((Tpk), i.e. the peak CO intensity at the line 
centre for all independent sightlines through the cloud, averaged over the projected cloud area) tend to be lower for young 
GMCs. There is no significant difference between the distributions of the total CO luminosity (Leo) or CO surface brightness 
(7co) however. The median Xco value for the young GMC sample {{Xco) = 6.9 x 10^" cm^^ (Kkms^^)~^) is ~ 50 per cent 
greater than the median Xco value for the entire MAGMA LMC cloud list. Some of this variation may be due to cloud size, 
since the discrepancy is reduced once we restrict our comparison to the control sample. The KS test results indicate that the 
difference is only marginally significant. 

For properties of the local interstellar environment, young GMCs appear to be distributed throughout the LMC in a similar 
fashion as star-forming GMCs (i.e. across a similar range of galactocentric radii and radial velocities), and they are detected 
across a comparable range of stellar mass surface densities and Hi column densities. There is a clear trend, however, for young 
GMCs to be located in regions where the Hi peak brightness ({rpk(HI))) is relatively low, and the radiation field (Go) is 
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relatively weak. In Section [7. II we argue that the latter reflects the fact that young massive stars are an important source of 
dust heating in LMC molecular clouds. 



5 CLOUD SCALING RELATIONS 

Empirical scaling relations between the basic physical properties of molecular clouds have become a standard test for poten- 
tial difl^erences between molecular cloud populations. Larson's initial work identified the three well-known "laws" obeyed by 
Galactic molecular clouds: i) a power-law relationship between the size of a cloud and its velocity dispersion, ii) a nearly linear 
correlation between the virial mass of a cloud and mass estimates based on molecular line tracers of H2 column density, which 



seemed to imply that molecular clouds are self-gra vitating and in app roximate virial balance, and iii) an inverse relationship 
between the size of a cloud and its average density ( Larson 1979I . 119811 ). S87 were subsequently able to measure the coefficients 



and exponents of the power-law relationships between the properties of 273 GMCs in the inner Milky Way, establishing the 
empirical expressions for Larson's Laws that have becom e the yardstick for studies of GMCs in other galaxies and in different 
interstellar environments (see e.g. BOS. lBlitz et al.ll2007r ). 



Most famously, the analysis by S87 showed that inner Milky Way clouds follow a size-linewidth relation of the form 



0.727?" 



km s , and that there is a strong, approximately linear correlation between the clouds' virial mass and their CO 



luminosity, Myir = 39Z/qo These two relations can be combined to provide expressions for the relationship between 

a cloud's luminosity and size, and luminosity and linewidth: Leo ~ 25-R^'^ Kkms"'^ pc~^ and Leo ~ ISOctv K kms~^ pc~'^. 
As noted by S87, clouds in gravitational equilibrium (M oc Ra^) with (Tv oc i?"'^ will follow a mass-size relation of the form 
M (X . This implies constant average mass surface density for molecular clouds, (EH2), which is related to the coefficient 
of the size-linewidth relation 
S87 found Co = 0.72 kms'^ 



Co = (Tv/a/R, via (Ehs) ~ 33IC0 Mq pc ^ for the cloud density profile that we have adopted. 
"■^ and (Shj) ~ 170 Mq pc~'^ for molecular clouds in the inner M ilky Way, sli g htly h igher 



pc 



Hever et all l|2008l ) have 



than the median surface density of the extragalactic GMCs analysed by B08 ((Ehj) ~ 130 M© pc" 
recently argued that the S87 estimate should be revised downwards to ~ 100 M0 pc~^, and we adopt this as our reference 
value in Figs. |4] to [8] 



5.1 The size-linewidth relation 

Fig. [4[a] presents a plot of velocity dispersion versus radius for GMCs in the MAGMA cloud list. The LMC clouds are clearly 
offset towards lower velocity dispersions compared to the R — relations derived by S87 and B08; some of the larger LMC 
clouds fall under the R — a^ relation determined by B08 (uv ~ 0.44^g ^gi?"'®*'''^) by a factor of ~ 3 in velocity dispersion. The 
discrepancy between the LMC clouds and the R~ relation derived by B08 should not be the result of cloud decomposition 
techniques, as we have identified cloud and parameterized CMC properties using the same algorithms as these authors. In 
the turbulent paradigm, the offset towards lower velocity dispersion at a given radius would suggest that the turbulent bulk 
motions within LMC molecular clouds are more quiescent than in the BOS GMCs. If, on the other hand, GMCs in the LMC 
manage to achieve rough dynamic equilibrium, then their relatively narrow linewidths would imply that they have lower mass 
surface densities than the B08 clouds. 



To fit the R — ay r e lation ship for the MAGMA data, we used the BCES bisector linear regression method presented by 
Akritas fc Bershadv ( 19961 ). which is designed to take measurement errors in both the dependent and independent variable, 
and the intrinsic scatter of a dataset into account^ For our analysis of all the Lars on relations, we use the bisector method 
because our goal is to estimate the intrinsic relation between the cloud properties (e.g. Babu fc Feigelsonll 19961 ). The best-fitting 
relation for all 125 clouds in the MAGMA cloud list, illustrated with grey shading in Fig. |3[a], is: 



logav = (-0.73 ± 0.08) + (0.74 ± 0.05) logR (BCES bisector). 



(4) 



Here, and for all other relations presented in this paper, the errors in the regression coefficients are derived using bootstrapping 
techniques; for relations determined using a BCES estimator, the y are c onsistent with the standard deviation of the regression 
coefficients derived according to eqn. 30 in lAkritas fc Bershadv l| 19961 ) unless otherwise noted. Although the errors in best- 
fitting relation that we derive using are small, we caution that the form of the relation depend s on the linear regression metho d 
(for a discussion of this issue in other astronomical contexts, see e.g. Tremaine et al. 2002 : Kelly 2007 : Blanc et al. 20091 ) . 
For comparison with other work, the best-fitting relations determined using the BCES ordinary least squares (henceforth 
OLS(Y|X)) and BCES orthogonal methods are: 



^ 'BCES' stands for bivariatc, correlated errors and intrinsic scatter. Software that implements this method is available from: 
[http : //www, astro .wise . edu/~mab/archive/stats/stats .html, For our analysis, we assume that measurement errors are uncorrelated. 
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Figure 4. Plots of (a) radius versus velocity dispersion; (b) radius versus virial mass; (c) radius versus CO luminosity; and (d) CO 
luminosity versus virial mass for the GMCs identified in the MAGMA LMC survey. In each panel, the light grey shaded area represents 
our BCES bisector fit to the 125 GMCs in the MAGMA cloud list and the l-cr uncertainty in this fit (see text). The black dotted 
line in each panel shows the standard relation for the S87 inner Milky Way data, while the black dashed line represents the relation 
determined for extragalactic GMCs by BOS. The dot-dashed grey lines represent constant values of CO surface brightness {Ico = li 
10, 100 Kkms~^, panel [b]); mass surface density (Ehj = 10, 100, 1000 Mq pc''^, panel [c]); and CO-to-H2 conversion factor {Xqo = 
0.4,4.0,40 X 10^'^ cm~2 (Kkms"^)^^, panel [d]). GMCs belonging to the complete MAGMA LMC cloud list are represented by small 
open squares, and GMCs in the high quality subsample are indicated using filled grey circles. The blue cross symbols represent GMCs 
without signs of active star formation. For comparison with the correlations presented in Fig. [5] to [S] Spearman's rank correlation 
coefficent and corresponding p-valuc for all clouds and the high quality subsample are indicated at the bottom right of each panel. 



log cTv = (-0.24 ± 0.09) + (0.40 ± 0.06) log R (BCES OLS(Y— X)) (5) 
logcTv = (-0.46 ± 0.12) + (0.56 ± 0.08) logi? (BCES orthogonal). (6) 
The best-fitting relation derived using the FITEXY estimator ( Press et al. 1992l l is: 



logcTv = (-0.37 ±0.17) + (0.50 ±0.11) log i? (FITEXY). (7) 

The intrinsic scatter of the MAGMA GMCs around the best-fitting i? — CTv relation is greater than the measurement errors 
in R and CTv across the observed range of cloud radii. The Co values of individual MAGMA GMCs vary between 0.21 and 
0.98 kms~^ pc~°'^. Although LMC clouds may be said to follow the same 7? — CTv relation as other extragalactic GMCs in 
the sense that the slope and amplitude of the derived best-fitting relations are similar, their molecular mass surface density 
is not strictly constant, but instead varies between 15 and 320 M© pc~^ (see Fig. [Hb]). An outstanding question, which 
we investigate further in Section [S] is whether the variation in is stochastic or whether it is related to changes in the 
environment of the GMCs. 
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5.2 The size-luminosity relation 

Fig. |4[c] presents a plot of CO luminosity versus radius for GMCs in the MAGMA cloud list. Although the Leo — R relation 
determined for extragalactic clouds by BOS (Leo ~ 7 .^'^^^^^R^-^^^^-^^) overlaps with the smaller MAGMA clouds, the slope 
of the Leo — R relation in the LMC is clearly shallower, such that the CO emission in large LMC molecular clouds is up to 
an order of magnitude fainter than for GMCs in nearby galaxies. A BCES bisector fit yields: 

log Leo = (1.39 ±0.12) + (1.88 ±0.08) log ii, (8) 

for the Leo — R relation. The relations determined using the BCES OLS(YjX), BCES orthogonal and FITEXY methods are 
also shallower than the B08 fit, with slopes of 1.67±0.09, 2.02±0.10 and 2.17±0.29 respectively 

The median CO surface brightness of GMCs in the MAGMA cloud list is 4.8 Kkms^^, and the median absolute deviation of 
the GMCs around this value is ~ 30 per cent. This corresponds to an average mass surface density of only (Eg,) ^ 20 M.p, pc"'^ 
for the MAGMA clouds if we adopt the Galactic Xeo value (2 x 10^° cm"^ (Kkms"^)"\ e.g. Istrong fc MattoxlHiiil ') and 



(Effa) ~ 50 Mq pc"^ if we use the median Xco value of the MAGMA cloud list (Xco = 4.7 x 10^° cm"^ (Kkms"^)"\ see 
Table[l}. Individual GMCs have Jeo values between 1.8 and 20.4 Kkms~^. Some of the variation in 7co appears to be related 
to the location of the GMCs: for clouds that are coincident with the stellar bar (defined as regions where E, > 100 M© pc~'^, 
see Fig.[T}, the median CO surface brightness is 9.1 K kms^^. We examine the relationship between the CO emission in LMC 
molecular clouds and the stellar mass surface density more closely in Section [6. 2. II 



5.3 The calibration between virial mass and CO luminosity 

Fig. |4[d] presents a plot of the virial mass estimate versus the CO luminosity for the MAGMA clouds. The LMC clouds 
appear to be in reasonable agreement with the slope of the relation determined for the B08 extragalactic CMC data (A/vir ~ 
7.6+2;6'f'co°*° °*)' ahhough offset slightly to higher Mvir values. A BCES bisector fit to the Mvir - Leo relation for the 
complete MAGMA LMC cloud list yields: 



log Mvir = (0.50 ±0.25) ±(1.13 ±0.06) log Leo. 



(9) 



The slopes of the relations determined using the BCES OLS(YiX), BCES orthogonal and FITEXY methods are 0.99±0.06, 
1.16±0.07 and 0.97±0.17 respectively. Even for the steepest BCES fit, the systematic variation in Xco with mass is small, 
corresponding to only a factor of two increase in Xco for CMC masses between 3 x lO" and 3 x lO'' Mq. The rms scatter of 
the MAGMA GMCs around the best-fitting relation corresponds to a similar variation in Xco- 



The median value of Xeo for the MAGMA LM C clouds is 4.7 x 10^" cm^^ (Kkms-i)-\ This is in excellent agreement with 
the values de rived by BOS and Blitz et al.l ([20071) for their extragalactic GMCs, and the LMC value obtained by the SEST Large 



Programme Jlsrael et al.ll2003l ). but lower than the value derived from the NANTEN LMC survey (7x10^'' cm-2 (Kkms 



Fukui et al.ll2008l ). While the CPROPS cloud decomposition algorithm aims to minimize instrumental effects, part of this 



discrepancy may be due to the difference in angular resolution between the two surveys. GMCs are constituted by dense 
(n « 10"^ c m~^) CO-bright peaks embedded in m ore diffuse gas with lower CO brightness; as noted by several previous 
authors (e.g. Pineda et al. 20091 : Bolatto et al. 20031 ) . observations with coarser resolution trace larger physical structures and 
hence derive larger values for Xco - In principle, the discrepancy could also arise from MAGMA's observational strategy, which 
excluded clouds in the NANTEN catalogue with low CO surface brightness (recall that Xeo = EH2/2.2/eo). In practice, 
however, we do not expect that our target selection has a significant impact on the average MAGMA Xeo value, since there 
are only four NANTEN clouds with Leo > 7000 Kkms"^ pc^ and peak Jeo < 1 Kkms"^ 



6 GMC PROPERTIES AND INTERSTELLAR CONDITIONS 

In this section, we explore the variation of the physical properties of the MAGMA LMC clouds in response to local interstellar 
conditions. We measure the local H i column density, stellar mass surface density, interstellar radiation field and external pres- 
sure for each GMC using the maps described in Section [2.21 taking the mean value of all independent pixels with integrated 
CO emission greater than 1 Kkms"'^. We measure the strength of correlations between the GMC and interstellar properties 
using the Spearman rank correlation coefficient, r, a non-parametric rank statistic that measures the strength of monotone 
association between two variables. The statistical significance of r is assessed by calculating the corresponding p- value, which 
is the two-sided significance level of r's deviation from zero. We consider p-values less than 0.01 to provide statistically sig- 
nificant evidence against the null hypothesis (i.e. that there is no underlying correlation between the variables). We regard 
|r| values greater than 0.6 as strong correlations (or anti-correlations if r < 0), \r\ values between 0.4 and 0.6 as moderate 
correlations, and \r\ values between 0.2 and 0.4 as weak correlations. We regard correlations with |r| values less than 0.2 to 
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Table 2. Results of the correlation tests between the properties of the MAGMA GMCs and properties of the interstellar environ- 
ment. The properties used in the comparison are listed in columns 1 and 2. The results of the correlation tests for the complete 
MAGMA cloud list are shown in columns 3 to 5; the results for the 57 high quality GMCs are shown in columns 6 to 8. The results 
list (r), (p) and ffp, where (r) is the median Spearman correlation coefficient obtained in the error trials (see text), (p) is the median 
of the corresponding p-values, and rjp is the standard deviation of the p-values. 
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be very weak and therefore unlikely to have practical significance, even if their p-value is small. 

Each correlation test was repeated for the complete MAGMA LMC cloud list and for the high quality subsample. Since the 
Spearman rank correlation test does not account for measurement uncertainties, we performed 500 trials of each correlation 
test in which we offset each cloud property measurement by a fraction of its uncertainty. As for the KS test in Section [H 
the value of a cloud property measurement in each trial was displaced by k/S.x, where Aa; is the absolute uncertainty in the 
cloud property measurement and is a uniformly distributed random number between -1 and 1. We consider our correlation 
result to be robust and significant if i) the correlation coefficients obtained in the trials are narrowly distributed around the 
original (i.e. unperturbed) value of r, ii) the corresponding p-values are narrowly distributed around zero, i.e. (p) ^ 0.01 and 
cjp ^ 0.01, and iii) \r\ ^ 0.2 for both the complete MAGMA cloud list and the high quality subsample. The results of all 
the correlation tests are presented in Table [2l As a reference for the values obtained in these comparisons, we conducted 
equivalent correlation tests for the Larson scaling relations shown in Fig. [31 and present the results in Table O We note that 
the size-linewidth relation is only a weak to moderate correlation: for all 125 GMCs in the MAGMA cloud list (r) = 0.48, 
while for the high quality GMCs (r) = 0.39. The size-luminosity and mass-luminosity relations are more strongly correlated 
((r) « 0.8); this is expected since Leo and Mvir are both dependent on R and a^. 



6.1 Comparison with Go 

In Fig. [5l we plot the radius, velocity dispersion, average peak CO brightness, CO surface brightness, mass surface density 
and CO-to-H[2 conversion factor of the MAGMA clouds as a function of the local interstellar radiation field. Go. All the plots 
show considerable scatter and none of the c orrelations satisfy our criteria for significance. Contrary to what might be expected 
from classic photodissociation models (e.g. Ivan Dishoeck fc Blacbil988i ). there is no general trend between the Go and Xco, 
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Table 3. Results of the correlation tests for Larson's scaling relations in the LMC. The relation is listed in column 1. The results 
of the correlation tests for the complete MAGMA cloud list are shown in column 2; the results for the 57 high quality GMCs are 
shown in column 3. The results are expressed in the form ((r), (p), (jp), where (r) is the median Spearman correlation coefficient 
obtained in the error trials (see text), (p) is the median of the corresponding p-values, and Cp is the standard deviation of the 
p- values. 
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even for clouds without signs of active star formation (i.e. where Go is dominated by the external field). An earlier analysis 
of the M AGMA data for the molecular ridge region also found that Xco was insensitive to variations in the radiation field 
strength (jPineda et al.ll2009l ). 



6.2 Comparison with the interstellar pressure 



In this subsection, we investigate whether the physical properties of the MAGMA clouds vary with the interstellar pressure. 
We estimate the total pressure at the molecular cloud boundary using the expression given by E89 for a two-component disc 
of gas and stars in hydrostatic equilibrium: 



(10) 



In this expression, ag and cr, are the velocity dispersions of the gas and stars respectively, Eg is the mass surface density of 
the gas and E* is the stellar mass surface density. The term in brackets on the right hand side of Equation [TU] is an estimate 
for the total dynamical mass within the disc gas layer. We assume a constant veloc ity dispersion of c tq = 9 kms~^ for the 



gas, based on the average dispersion of the Hi line profiles across the LMC (see also IWong et al. 200S), and a* = 20 kms 



van der Marel et al.ll2002l '). For E*, we use the stellar mass surface density map of I Yang et al.l (|2007l ') described 



for the stars 

above, and we estimate Eg 

conversion between A''(Hi) and Eg includes a factor of 1.36 by mass for the presence of helium). We note that the Hi and 
stars make similar contributions to the total mass surface density at the locations of GMCs in the LMC: (E«) = 48 Mq pc~^, 



Ehi[M0Pc" ] = 1.089 X 10" °iV(HI)[cm"^] directly from the Hi column density map (the 



and (Eg) = 30 Mq pc" 
the Milky Way disc. 



for comparison, E89 adopted (E,) = 55 Mq pc"^ and (Eg) = 12 Mg pc as characteristic values in 



Along sightlines to GMCs, our definition of Ph provides an estimate of the pressure at the surface of the molecular cloud 
rather than at the disc midplane, since the weight of a GMC is likely to make a significant contribution to the total midplane 
pressure. Further caveats are (i) the assumption that the disc is supported against gravity solely by the observed velocity 
dispersions, and (ii) that the interaction between the LMC and the Galactic halo may make a significant contribution to 
the pressure that is not accounted for by our simple estimate in Equation 1101 In the following, we first compare the cloud 
properties to two independent components of Ph for which we have an empirical tracer: i) the stellar mass surface density, E* 
and ii) the atomic gas surface density. Eg. 



6.2.1 Comparison with E, 



In Fig. [6l we plot the properties of the MAGMA clouds versus E*. The correlation tests indicate weak but significant 
correlations between E* and (Tpk) (panel [c]) and Ico (panel [d]); these correlations are strengthened if we restrict our 
analysis to high quality GMCs. The existence of these trends suggests that the GMCs in the LMC do respond to the presence 
of the galaxy's stel l ar population, despi te evidence that the stellar bar may be physically offset from the gas disc (e.g. 
Zhao fc EvaniboOOl : iNikolaev et aklbood ). 



6.2.2 Comparison with Eg 

In Fig. [T] we plot the cloud properties as a function of the atomic gas surface density. Eg , which is estimated from the total 
Hi column density along the line-of-sight. We find a weak but robust correlation between Eg and ctv (panel [b]); again, the 
correlation is stronger if only high quality GMCs are considered. There is some indication that Ehj (panel [e]) and Xco 
(panel [f]) also increase with the H i column density - as would be expected if o"v increases without a corresponding increase in 
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Figure 5. Properties of the MAGMA GMCs, plotted as a function of the local interstellar radiation field, Go/Go,©: (a) radius; (b) 
velocity dispersion; (c) average peak CO temperature; (d) CO surface brightness, Ico\ (s) molecular mass surface density, SH2i ^.nd (f) 
CO-to-H2 conversion factor, Xqo- The plot symbols are the same as in Fig.|4l The horizontal dotted lines in panels [e] and [f] indicate 
= 100 M0 pc-2 and Xqo ~ 2 x 10^° cm'^ (Kkms~l)~l, i.e. values that apply to GMCs in the inner Milky Way. The Spearman 
rank correlation coeffient and corresponding p-value for the complete MAGMA cloud list and high quality subsample are indicated at 
the top left of each panel. 



cloud size or CO luminosity - but the correlations are not significant if only high quality GMCs are considered. An important 
caveat for interpreting Fig. [7] is that Hi line profiles in the LMC are notorio usly complex, exhibit ing two well-defined velocity 
components in some parts of the LMC's disc (especially in the south-east, e.g. Luks &: Rohl3l992l '). The CO emission is almost 



invariably associated with only one Hi velocity component (|Wong et al.l 120091 '). so the Hi column density that is physically 



associated with a CMC is almost certainly overestimated by the total A'^(Hi) along the line-of-sight in these regions. Excluding 
GMCs located in the south-east of the LMC (i.e. clouds with right ascension above 05h38m and declination below -68d30m 
(J2000)) from our comparisons between the cloud properties and Eg suggests that the correlations between Eg and (t„ is not 
severely contaminated by the contribution of the secondary Hi velocity component to Ai'(Hi). 
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Figure 6. Properties of the MAGMA clouds, plotted as a function of the stellar mass surface density, S». The panels, plot symbols and 
annotations are the same as in Fig. [5] The grey shading indicates regions where S* ^ 100 Mqpc"'^, corresponding to the stellar bar 
region denoted by the ellipse in Fig. [T] 



6.2.3 Comparison with Ph 

The comparison between properties of the MAGMA clouds and Eg is helpful for interpreting the plots in Fig. [S] In particular, 
both (Tv and Ehj increase in regions with higher P^, following their behaviour in regions with high Hi column density. The 
dominant physics underlying these correlations may therefore be the importance of an atomic shielding layer for the survival 
of H2 molecules, rather than pressure regulation. and Ph is more strongly correlated with Ph than with Eg however, and 
we observe that there are few clouds with low Ico or Eh2 in regions with high E, (panels [d] and [e] of Fig. [6l). Insofar as Ico 
and Eh2 are reliable tracers of the H2 surface density, this provides some indication that shielding alone may not regulate the 
H2 surface density in regions of the LMC where the interstellar gas pressure is high. 



© 0000 RAS, MNRAS 000, 000-000 



18 A. Hughes et al. 



2.5 



O 

£j 2.0 
"m 

1.5 

o 



I ' ' ' I ' ' ' I ' ' ' I 

1^ [a] rail = (0.12,0.19) 
rhq = (0.28,0.04) 



• □ 8 



1.0 



0.4 



0.3 



0.2 




0.1 

V 



0.0 

-0.1 

-0.2 
3.0 



I ■ ' ' I ■ ' ■ I ' ' ■ I 
[c] Tall = (-0.07,0.36) 
rhq = (-0.15,0.24) 




I 2.5 
O 



2.0 



O 



1.0 



I ■ ' ■ I ■ ' ■ I ■ ' ' I ' ' ' I 

[e] raii = (0.25,<0.01) 
rhq = (0.26,0.06) 




I . . I 



1.0 1.2 1.4 1.6 1.8 

iog(i:g/[Msoi pc-2]) 



1.0 1.2 1.4 1.6 1.8 2.0 



1.0 1.2 1.4 1.6 1.8 2.0 



2.0 



1.0 
0.8 
0.6 - 



I 



> 

'bC 0.2 
0.0 
-0.2 

_ 1.5 



M 1-0 



I ■ ' ■ I ■ ' ■ I ■ ' ' I ' ' ' I 

[b] Tall = (0.30,<0.01) 
rhq = (0.43, <0.01) 




I . . I . I . . I 



1.0 1.2 1.4 1.1 



1.8 2.0 



o 
u 

o 



0.5 



0.0 



I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 
[d] ran = (-0.05,0.58) 
rhq = (-0.03,0.90) 




_L 



I ■ ■ I ■ I ■ ■ I 



1.0 1.2 1.4 1.6 1.1 



I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 
m rail = (0.28,<0.01) 
rhq = (0.29,0.03) 



log(Eg/[Msoipc-2]) 



2.0 




Figure 7. Properties of the MAGMA clouds, plotted as a function of the atomic gas surface density, Eg. The panels, plot symbols and 
annotations are the same as in Fig. \5\ 



7 DISCUSSION 

7.1 Properties of young GMCs 

In Section m we investigated the physical properties of the young GMC candidates in the MAGMA LMC cloud list, finding 
that sightlines through these clouds tend to have lower Tmax and (Tpk) than sightlines through star-forming GMCs of equiv- 
alent size. Models of the^^CO{J = 1 0) emission f rom GMCs in dicate that the emission arises from a large number of 



optically thick clumps that are not self-shadowing (e.g. IWolfire et al.lil99 3). For each independent sightline through a GMC, 



the observed peak CO brightness is then a measure of the total projected area of the optically thick clumps within the beam 
area, weighted by their brightness temperatu re, which should c orrespond to the kinetic temperature of the CO-emitting gas if 
emission in the clumps is optically thick (e.g Malonev &: Blac^ 1988, ). The usual assumption is that the distribution of clump 



sizes and the brightness temperature do not vary significantly between sightlines, and that the observed CO brightness mea- 
sures the number of clumps - and hence the total amount of molecular gas - within the telescope beam. While this assumption 
may be justified for clouds in the inner disc of the Milky Way (S87), it is worth noting that the peak CO brightness results 
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Figure 8. Properties of the MAGMA clouds, plotted as a function of the interstellar pressure, P^. The panels, plot symbols and 
annotations are the same as in Fig. [5] 



from a combination of the brightness temperature, the number of CO-emitting clumps within the beam, and their average 
size. A possible interpretation for the different average peak CO brightness of star-forming and non-star-forming GMCs is 
that there are fewer CO-emitting clumps in GMCs without star formation, leading to a lower angular filling factor of CO 
emission. Alternatively, colder gas temperatures in the dormant CO-emitting substructures within non-star-forming GMCs 
could lead to a lower average brightness temperature for the clump ensemble. 

As our young CMC candidates tend to have a lower (Tpk) and Tmax values than star-forming GMCs of comparable size, it 
might be expected that their total CO luminosity would also be fainter. However, the KS tests reveal no systematic differ- 
ences between the distributions of Leo and Ico for the various cloud samples, which suggests that the regions of high CO 
brightness in the star-forming GMCs are restricted to a relatively small number of pixels in the MAGMA data subcubes. 
This is consistent with the view of star formation as a highly localised process: occupying only a small fraction of the total 
cloud volume, star-forming clumps have temperatures and densities that are much higher than in the bulk of the CMC, most 
of which does not participate directly in star formation. The lower (Tpk) and Tmax of non-star-forming GMCs would seem 
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Figure 9. The radiation field at the location of the MAGMA clouds, plotted as a function of (a) the Ho surface brightness, as observed 
by SHASSA iGaustad et al.: 2001), and (b) the stellar mass surface density. Plot symbols are the same as in Fig. |4] The Spearman 
correlation coefficient and corresponding p- value for all 125 MAGMA GMCs and the high quality subsample are indicated at the top loft 
of each panel. 



to provide some preliminary evidence that the general characteristics of substructure within non- star- forming GMCs in the 



Milky Way - i.e. coole r clumps that are less massive and more diffuse than in star-forming GMCs (e.g. IWilliams et al.lll994 
Williams fc Blitzlll99 - wiU also be found to apply in the LMC. 



Finally, we found a clear trend for the young CMC candidates to be located in regions where Go is relatively weak. A 
straightforward explanation for this result is that once young massive stars begin to form, they are an important source of 
dust-heating within LMC molecular clouds. Importantly, this would mean that Go is not strictly tracing the ambient (i.e. 
external) radiation field incident on a CMC, but instead contains a signficant contribution from the massive stars that are 
within - or have recently emerged from - their CMC progenitor. A strong correlation between Go and Ha surface brightness 
for all the MAGMA LMC clouds (see Fig. [5| indicates a strong connection between dust heating and star-formation activity, 
which would seem to support this interpretation. 



7.2 The origin of Larson's laws 



Despite the longevity of Larson's scaling relations, a complete theoretical explanation for the origin of the size-linewidth 
relation is still lacking. Considerable effort has been devoted to de monstrating that the observed size-linewith relationship can 



be reproduced by realistic models of interstellar turbulence (e.g. iMac Low fc Klessenll2004l : lBallesteros-Pare des_et^ 



2007 



and references therein), and that turbulence in the cold gas phase is both universal and nearly invariant fe.g. iHever fc Brunt 
200J), but differences between the working variables and descriptive tools available to theory and observation represent a 
significant obstacle to making turbulent models empirically falsifiable. In many instances, however, older explanations for the 
origin of Larson's laws have not been thoroughly tested, since there are few extragalactic datasets that are comparable to the 
Milky Way molecular cloud samples that guided the development of these early models. In this section, we therefore consider 
whether the MAGMA data are consistent with two analytic models for the origin of Larson's laws: the photoionization- 
regulated star formation theory proposed by M89 and the model of molecular clouds as virialized polytropes proposed by 
E89. 



7.2.1 Comparison with\McEm uMm) 



The photoionization-regulated theory of star formation proposed by M89 provides one possible explanation for the origin 
of Larson's laws. In this theory, molecular clouds evolve towards an equilibrium state where energy injection from newborn 
low-mass stars halts the clouds' gravitational contraction. Equilibrium depends on the level of photoionization by the inter- 
stellar far-ultraviolet (FUV) radiation field because ambipolar diffusion governs the rate of low-mass star formation in the 
cloud; in turn, the rate of ambipolar diffusion is regulated by the ionization fraction, which depends on the interstellar FUV 
radiation field in the bulk of the molecular cloud. A second factor that determines the equilibrium state of the GMCs is 
the local dust abundance, since dust shields H2 molecules against photodissociation. Notably, the M89 theory predicts that 
molecular clouds in equilibrium should have uniform extinction, rather than constant column density: clouds in environments 
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with low dust-to-gas ratios and/or strong radiation fields will require larger column densities to attain the equilibrium level 
of extinction, and should therefore follow a, R — relation with a higher coefficient than the Milky Way relation. Recent 
measurements of the surface densities of extragalactic GMCs contradict this prediction, instead showing that molecular cloud 
surface densities in lo w metallicity envir onments are similar to, or even lower than, the average mass surface density of Milky 
Way clouds (e.g. BOS, iLerov et al.ll2007l ). The median mass surface density of the MAGMA clouds is ~ 50 M0 pc^^, in line 
with these results. 



It is worth noting, however, that low mass surface densities (i.e. less than ~ 100 M© pc~^ for GMCs in sub-solar metallicity 
environments) are not necessarily inconsistent with the M89 model: the complete prediction by M89 is that the velocity dis- 
persion of molecular clouds should increase as the dust-to-gas ratio decreases, provided that the densities of the CO-emitting 
clumps within GMCs are comparable to the densities of the clumps in Galactic molecular clouds, and that appropriate cor- 
rections for the angular filling fraction of the CO-emitting clumps, C, have been applied. If the CO-emitting clumps within 
GMCs have C ~ 1, densities of nn ~ 10'^ cm~'^, and clump-to-cloud extinction ratios of Ac/A^ ~ 0.3 - values that M89 
regards as typical for molecular clouds in the inner Milky Way disc - then the photoionization-regulated star formation theory 
predicts that equilibrium is achieved for visual extinctions between A„ ~ 4 and 8 mag. Equations 5.5 to 5.7 in M89 show 
that an equilibrium extin ction value oi Ay ~ 1 can none the less be obtained in an environment with the LMC's typical 
dust-to-gas ratio (~ 0.3, Dobashi et al. 20081 ) if the CO-emitting clumps within GMCs have an angular filling fraction of 
C — 0.25 and a typical density of nu = 10*'^ cm~^, and if the fraction of the total cloud mass residing in these dense clumps 
is ~ 45 per cent. While agreement with the M89 model is theoretically possible, it requires the density contrast in GMCs in 
the LMC to be more extreme than in Milky Way clouds. C values less than unity have previously been invoked to explain the 
low Tpk measurements for LMC clouds (e.g. IWolfire et al.lll993 : Kutner et al.lll997l : ICarav et al.ll2002l ). but average densities 
of nu ~ 10^'^ cm~^ for the CO-emitting clumps seem less plausible. Excitation analyses of millimetre and submillimetre 
spectral line observations (e.g. Heikkila et al. 19991 : Pineda et al. 20081 : Minamidani et al. 20081 ) have reported clump densities 
between 10* and 10® cm~^, but as these studies targeted the LMC's brightest star-forming regions it remains uncertain 
whether similar clump densities would be common in molecular clouds throughout the LMC. 



A further problem for the photoionization-regulated star formation theory is that there is no sign of a correlation between Go 
and or Ico for GMCs without signs of star formation (blue crosses in panels [d] and [e] of Fig. [5]). In the scenario postu- 
lated by M89, equilibrium is only achieved prior to the onset of massive star formation: after this, young massive stars rapidly 
disrupt their natal clouds, destroying the relationship between the gas column density, dust abundance and ambient radiation 



field. The GMCs that are designated as star-forming in the MAGMA sample contain at least one O star (jKawamura et al. 



2OO9I ) and their young stellar content makes a significant contribution to the radiation field within the cloud (see Table [T]). For 
GMCs without star formation, however, and Ico should show a correlation with Go if their column density is regulated 
by the ambient radiation field, assuming that the dust-to-gas ratio is roughly constant in the environments of non-star-forming 
GMCs across the LMC. No such correlation is apparent in panel [d] or [e] of Fig. [5] 



Finally, we note that a correlation between the internally-generated radiati on field and the H 2 mass surface density would 
be expected if the star formation efficiency of molecular gas were constant ( Lerov et al .' 2008), since GMCs with higher H2 



column densities - and presumably higher volume densities - should have a higher surface density of star formation. The 
correlation tests in Section [6] indicated that E^r^ for star- forming GMCs in the high quality subsample are associated with 
higher values of Go, but the trend is not significant if we consider all the GMCs in the MAGMA cloud list, or if we use Ico 
rather than Eh2 to trace the H2 mass surface density. 



7. 2. 2 Comparison with Elmeareen , 1989 ) 



Another explanation for the origin of Larson's laws was put forward by E89, who proposed that molecular clouds and their 
atomic envelopes are virialized, magnetic polytropes (i.e. with internal pressure P that varies with the density p according 
to P = Kp", where A' is a constant and n is the polytropic index) and with external pressure that is determined by the 
kinetic pressure of the interstellar medium. In this theory, the radius and density of a molecular cloud complex adjust so 
that the pressure at the boundary of the atomic envelope equals the ambient kinetic pressure. In this case, molecular clouds 
within a galaxy have a similar mass surface density because the ambient pressure throughout the galaxy is also roughly 
unif orm. While the notion that molecular clouds can achieve dynamical equilibrium within their lifetime has been disputed 
(e.g. Hartmann et al.l 200 ll). recent work has highlighted the potential importance of pressure for the formation of molecular 



gas (e.g. IWong fc Blitzlbood : iBlitz fc Rosolowskvl[2cio3 . bood ) 



The molecular cloud model proposed by E89 predicts that Co and Ehj should scale with the total mass of the atomic-|- molecular 
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Table 4. A summary of the model parameters which are substituted into Equation [TJ] in or der to produce t he plots in Fig. IIUI The first 
column describes the metallicity: either a shallow radial gradient of —0.05 dex kpc"'^ (e.g. ICole et al.ll20oj ). or a constant Z/Zq = 0.5 
for all clouds in as sumed. We adopt 05hl9m30s -68d53m (J2000) for the location of the LMC's kinematic centre, and i = 35° for the 
LMC's inclination l lWone et al.ll2009h . The second column lists k, the factor by which the measured value of Gq/Gq^q has been reduced 
for star-forming clouds in the MAGMA cloud list. We introduce k to account for the fact the radiation field in Equation 1121 is the 
external field incident on the molecular cloud, i.e. it does not include the contribution from young massive stars within the cloud. The 
third column describes the relationship between the total atomic-|- molecular mass, Mt, and the molecular mass, Mm = M^\-^, that we 
adopt for the clouds. For Ml to M4, we use the mass dependence adopted by E89. The mass dependencies for M5 and M6 are constructed 
such that the total mass of the cloud complex increases more rapidly than the molecular mass; the coefficient of the mass dependence 
for these models is then chosen such that Mt ^ Mm for the observed range of MAGMA GMC masses. 



Model Identifier 


Z/Zq 


k 


Mass Dependency 


Ml 


0.5 


1.0 


Mt = 2Mm 


M2 


0.5 


2.0 


Mt = 2Mm 


M3 


Gradient 


1.0 


Mt = 2Mm 


M4 


Gradient 


2.0 


Mt = 2Mm 


M5 


Gradient 


1.0 


Mt/[MQ] = {Mm/21A5 Mq)1-5 


M6 


Gradient 


1.0 


Mt/[MQ] = (M„/141.42 M0)2-O 



cloud complex, due to the increasing weight of the atomic gas layer surrounding large molecular clouds. In Sections 15.11 and [6l 
we saw that the scatter in the Larson-type scaling relations for the MAGMA clouds implies order of magnitude variations 
in Co and Eh2, and also evidence for a weak correlation between and the interstellar pressure (Fig. [8le]). These trends 
motivate us to examine more closely whether the MAGMA clouds are consistent with the E89 theory for the origin of Larson's 
Laws, once variations in the local radiation field, metallicity and mass of atomic gas surrounding the GMCs are taken into 
account. 

In lieu of Larson's third law, E89 predicts that the total (i.e. atomic-|~molecular) mass surface density of molecular cloud 
complexes will vary according to 

A/f / P \ "^/^ 

-I ^ 190 ± 90 Mopc-^ (11) 



i?2 Vl04kBcm-3K/ 

where Pe is the external pressure on a molecular -|-atomic cloud complex, Alt is the total mass of the atomic-|-molecular 
complex and Rt is the radius of the complex. This equation can be written in terms of observed molecular cloud properties: 

Pi:"''°^''H l04kBcm-3K j [t¥m^) ""^P' ' ^''^ 

where Rm and Mm are the radius and mass of the molecular part of the cloud complex. The range of coefficients in these 
equations arise from solutions to the virial theorem for plausible values of the adiabatic index and variations in the ratio of the 
magnetic field pressure to the kinetic pressure (see fig. 3 in E89). To explore the relative importance of the parameters that 
contribute to the right hand side of Equation 1 121 we defined six basic models that use different estimates for the total mass, 
metallicity and external radiation field of the molecular cloud complexes. For all models, Pe is assumed to be the same as the 
kinetic pressure of the ambient interstellar medium, i.e. Pe = Ph/(1 + a + /3) where q ~ 0.4 and /3 ~ 0.25 are the relative 
contribution by cosmic-rays and the magnetic field to the total pressure, and Ph is estimated according to Equation 1 101 The 
model parameters are summarized in Tabled) 

In Fig. llOl we plot the observed values of Ehj for the MAGMA clouds against the values predicted by Equation [12] for each of 
the six models. Fig. [10] shows that the measured values for Ehj for the MAGMA clouds are significantly lower than predicted 
by E89, regardless of which model we adopt. The average mass surface density is sensitive to variations in the metallicity 
and radiation field, but the LMC would need to have Z Zq and Go ^ Go,© for the MAGMA clouds to be consistent with 
the closest line of equality in Fig. [TD] (which corresponds to the lowest value of the coefficient in Equation I12|) . The models 
in which the total complex mass increases faster than the molecular mass (M5 and M6) demonstrate a better agreement 
with the slopes predicted by E89, but these models are not physically realistic since they imply that the mass of the Hi 
envelopes surr ounding the MAGMA clou ds is almost an order of magnitude greater than the total Hi mass of the LMC 
(5 X 10** Mq, Istavelev-Smith et al.|[20o3 ). Reducing the coefficient of the mass dependency in these models would lower the 



total H I envelope mass and also shift the data towards the line of equality in Fig. 1121 but it leads to the unphysical solution 
that Mt < Mm for low mass MAGMA clouds. 

The total atomic+molecular mass of the cloud complex is perhaps the most important source of uncertainty in the present 
comparison, and we are currently undertaking an analysis of the H i and MAGMA CO datasets that should place empirical 
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Figure 10. The observed values of for the MAGMA clouds versus the values predicted by E89 for the models of the metallicity, 
external radiation field and total cloud mass presented in Table |4] Plot symbols are the same as in Fig. (4] The solid line represents 
equality between the observed and predicted values for the coefficients adopted by E89. The dashed lines indicate the potential shift in 
the line of equality for the maximum range of the coefficients in Equation 1121 



constraints on the mass of the Hi gas that is associated with individual MAGMA clouds. However, two further comments 
regarding the Hi data and its potential to constrain the E89 model are worth noting here. First, visual inspection of the CO 
and Hi LMC maps shows that the majority of GMCs are not isolated, but are instead located within a group of molecular 
clouds that appear to share a common Hi envelope. This possibility is not explicitly addressed by E89, but it is important 
insofar as neighbouring molecular clouds will partially shield each other, reducing the total amount of Hi required for the 
shielding layer. Second, we note that CO emission is detected in region s where the total Hi column density through the LMC 
is only Ai'(HI) « 1 — 2 x 10 cm" -2 (refer fig. 2 of I Wong et Zll2009[ ). This is problematic insofar as the H2 self-shielding 
theory invoked by E89 requires a uniform shielding layer with A'^(HI) = 1.2 x 10^^ cm~^ for a molecular cloud situated in an 
ambient radiation field Go/Go, o = L5 and metallicity Z/Zq — 0.5, assuming an density of n « 60 H cm~'^ in the atomic layer 
( Federman et al. 197^). A'^(HI) = 1.2 x 10^^ cm~^ is probably a lower limit: if the gas in the shielding layer is clumpy, then the 
required column density increases by a factor ///q, where / — {n)^/{n^) is the volume filling factor. Part of this discrepancy 
may be because sightlines towards GMCs in the LMC have significant optical depth: Hi absorption studies indicate peak 
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optical depths between r = 0.4 and 2.0 near regions with C O emission, implying that the true Hi column density is a factor 
of t/[1 — e(—r)] — 1.2 — 2.3 greater than the observed value ( Dickey et al. 1994 : Marx-Zimmer et al. 2000l ). A recent study of 
the LMC's far-infrared emission has also proposed a widespr ead cold atomic gas c omponent with significant optical depth in 
the LMC in order to explain an excess of emission at 70 /xm ( Bernard et al. 2008h . 



More generally, we note that plausible variations in the interstellar pressure due to e.g. a widespread ionized gas component 
or a high rate of cosmic-ray escape would not have a significant impact on the plots in Fig. 1101 since the exponent of the 
pressure term in Equation [12] is small. A significant mass of H2 without CO emission associated with each CMC would tend to 
rather magnify than reduce the discrepancy between the MAGMA clouds and the E89 model, moreover, except in a physically 
improbable scenario where the CO-dark H2 gas has a greater average density than the CO-emitting region of the GMC. One 
possible resolution may be the porosity of the LMC's interstellar medium: as noted above, the calculations in E89 assume a 
uniform gas layer which provides more effective shielding than gas that is highly clumped. A comparison between the volume 
filling factor of H i in the solar neighbourhood and in the LMC, plus the inclusion of clumpy cloud structure in the E89 theory, 
would be required to test this hypothesis however. 



7.3 GMC properties: trends with environment 



Previous comparative studies of the Milky Way and nearby extragalactic GMC populations have found that GMC properties 
are relatively uniform, and mostly insensitive to variations in env ironmental parameters such as the radiation field and dust 



abundance (e.g. BOS, Rosolowskv et al. 20031 : Rosolowskv 2007). A possible explanation for this result is that GMCs are 



strongly bound and hence largely decoupled from conditions in the local ISM. BOS have cautioned, however, that the apparent 
universality of GMC properties as measured from CO observations may simply reflect the physical conditions required for 
CO emission to be excited, and that the properties and extent of the H2 material surrounding a GMCs high density CO- 
emitting peaks might indeed be sensitive to local environmental factors, a conc l usion that is supported by far-infrared studies 
of molecular clouds in nearby dwarf galaxies (e.g. Israeli 1997 : Lerov et al. 2007 . 20091 ) . Our analysis tends to support the view 
that CO-derived properties of GMCs are mostly insensitive to environmental conditions, but there are several exceptions that 
we discuss below. 



7.3.1 GMC properties and T.:t 



In Section [6.2.11 we showed that (Tpk) and Ico increase for GMCs located in regions of high E», without a corresponding 
variation in i? or ctv. A key variable to explain this result may be the role of the old stars in dust production, as a recent 
analysis of Spitzer mid- and far-infrared data for the LMC indicates tha t VSGs and polycyclic aromatic hydrocarbon (PAH) 
molecules are overabundant in the LMC's stellar bar (Parad is et al. 200^). More specifically, if the abundance of VSGs and/or 
PAHs were increased by the ejecta from mass-losing old stars, then photoelectric heating of the molecular gas may be more 
efficient in regions with high E,, potentially raising the CO excitation temperature. Alternatively, a higher overall dust abun- 
dance could reduce the selective photodissociation of CO molecules and lead to a higher abundance of CO relative to H2 for 
GMCs in the stellar bar. This could occur if H2 readily sel f-shields, while the surv ival of CO molecules relies more on the 
attenuation of the photodissociating radiation by dust (e.g. Malonev fc Blac^ igssl ). Although the CO emission from Milky 
Way molecular clouds is approximately independent of variations in the CO abundance - firstly due to the high optical depth 
of the ^^CO(J = 1—^0) line, and secondly b ecause the angular filling factor of the CO-emitting clumps at any particular 
radial velocity within t he cloud is ~ 1 (see e.g. IWolfire et al]|l993l ) - this may not be true of LMC molecular clouds (see e.g. 
Malonev fc Blacklll988l ). The low values of Tpk for i^CO(J = 1^0) emission in LMC clo uds suggests that the a ngular filling 
facto r of the CO-emitting clumps is indeed relatively small (see e.g. Kutner et al. 1997; Johansson et al. 19981 : Garav et al 
20021 ). in which case higher CO-to-H2 ratios might produce higher values of (Tpk) and 7co. 



Contrary to what might be expected if higher interstellar pressures promote the formation of molecular gas, there is no 
evidence that R or ctv for the LMC molecular clouds increases in regions with high E* . The absence of these correlations is 
perhaps not very significant, however, since our estimate for Ph is dominated by Eg for regions with E. sC 60 M0 pc"^; if the 
Hi emission along sightlines towards GMCs has significant optical depth. Eg will dominate to even higher E, thresholds 
(E, ~ 100 M0 pc~^ for r = 1). A correlation between Go and E* might be expected if regions with str ong stellar gravit: 



promoted the formation of overdensities within GMCs and enhanced the local star formation rate (e.g. [Lerov et al.[ [20081 ) 
or, alternatively, if old stars made a significant contribution to the dust heating. Stronger external heating in the stellar bar 
region has previously b een invoked as an explanation for the higher ^^CO(J = 2 — >■ l)/( J = 1 — >■ 0) intensity ratios in inner 
LMC molecular clouds ([Sorai et al ] |200ll ). While the correlation tests in Section[6]provided no clear evidence for a relationship 
between Go and E*, we do not observe clouds with low Go values located in regions with E* ^ 100 M© pc~^ (see Fig.[9[b]), so 
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it is possible that enhanced star formation due to strong steUar gravity and/or dust heating by old stars becomes important 
at higher stellar densities. 



7.3.2 GMC properties and N (H i) 

In Section r6.2.2l we found a weak but significant correlation between A'^(Hi) and ctv (Fig- Ob])- Its interpretation is somewhat 
uncertain, however, since there is no obvious trend between A''(Hi) and the GMC radius, which would be expected if larger Hi 
column densities were associated with more massive GMCs (assuming that the average density of GMCs also remains constant 
across the LMC). Some positive association between A^(IIi) and R would be expected, moreover, simply as a consequence of 
the size-linewidth relation for LMC molecular clouds. Perhaps the simplest explanation is that the mass and size of GMCs 
increase with Ai'(Hi), but that the GMC radius derived from CO observations is not a sensitive tracer of the true cloud 
size. Another possibility is that the density of molecular clouds genuinely increases with the local H i column density; this 
could account for some of the scatter in the R — relation but would imply tha t Sht and Xco - assuming the averag e 



CO brightness temperature is constant - should also increase with A^(IIi) (see e.g. Dickman et al.l 19861 : Hever et al 



2001 



A third possibility is that the degree of virialisation in GMCs varies across the LMC, i.e. GMCs in regions with high A'^(H i) 
are less gravitationally bound than GMCs in regions with low A''(Hi), due to the higher external pressure. Provided that 
the average CO brightness temperature of GMCs in the LMC remains constant, we would also expect the observed value of 
Xco to increase with A'^(IIi) in this case. Our analysis in Section [6] revealed some evidence for trends between Ehj, Xco and 
A''(Hi), although the correlations involving the high quality subsample did not satisfy our criteria for significance. There is 
sufficient scatter in Figures |4] and [7] that none of these explanations can be definitively ruled out. 



7.3.3 GMC properties and Ph 

Although the MAGMA clouds do not follow the predictions of the E89 molecular cloud model, we find that CTv and Ehj 
increase with the interstellar pressure (panels [b] and [e] of Fig. [8]). These correlations suggest that the external pressure on 
a GMC in the LMC may indeed play a key role in regulating its dynamical properties, although we caution that we cannot 
readily distinguish between the role of pressure and shielding in our analysis due to the dominant contribution of Eg to our 
estimate for Ph. The low absolute values of Ehj and Ico for the MAGMA clouds imply that their average internal pressure 
is als o low. From the virial theo rem, a GMC with ~ 50 M© pc^^ will have internal pressure Ph/ks ~ 5 x 10* K cm~^ 
(e.g. iKrumholz fc McKeelliooil ) : this is not much greater than the average external kinetic pressure for the MAGMA GMCs, 



{Ph/ks) ~ 3.9 X 10* K cm~^, estimated from Equation 1101 Plausible optical depth corrections for the Hi emission along 
sightlines towards GMCs (i.e. for t values between 0.4 and 2) would effectively balance these estimates for the average pressure 
internal and external to the cloud, although the dense clumpy structure within the GMC would be likely to remain significantly 
overpressured. A gentle pressure gradient across the molecular cloud boundary may explain why we have found that some 
properties of the MAGMA clouds are correlated with environmental conditions. It would be interesting to test whether GMCs 
in nearby spiral galaxies - which have (Sh^) ~ 1 30 M© pc~^ (B08) and presumably higher internal pressures relative to the 
surrounding ISM (see also Krumholz et al. 2009l l - exhibit any of the correlations that we have identified for the MAGMA 



GMCs. 



8 CONCLUSIONS 

This paper has presented Mopra Telescope observations of ^^CO(J = 1 — >■ 0) emission from a sample of 125 GMCs in 
the Large Magellanic Cloud. The data described here were obtained as part of MAGMA, an ongoing mapping survey of 
the '^CO(J =1-^0) emission from molecular gas in the Magellanic Clouds. The MAGMA data, which have an angular 
resolution of 45 arcsec and a spectral resolution of O.lkms"^, will be made available to the astronomical community upon 
completion of the survey, and should be a rich resource for follow-up studies with millimetre and submillimetre facilities on 
the Atacama Plateau. A discussion of MAGMA's observational strategy, data products and a final GMC catalogue will be 
presented in a future paper (Wong et al., in preparation). In this paper, we examined the empirical scaling relations between 
the basic physical properties of the GMCs in our current LMC cloud list, and dependencies between the cloud properties and 
the local interstellar environment. We report the following results and conclusions: 

1. The observed GMCs have radii ranging between 13 and 160 pc, velocity dispersions between 1.0 and 6.1 kms~^, peak 
CO brightnesses between 1.2 and 7.1 K, CO luminosities between 10"^'^ and 10^'^ K kms~^ pc^, and virial masses be- 
tween 10*'^ and 10^ * M©. The clouds tend to be elongated, with a median major-to-minor axis ratio of 1.7. These values 
are comparable to the measured properties of Galactic GMCs. The average mass surface density of the observed clouds is 
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50 M0 pc , approximately half the value determined for GMCs in the inner Milky Way catalogue of lSolomon et al.l (119871 1 



2. The MAGMA cl ouds exhibit scalinR relations that are sim ilar to those previously determined for Galactic and extragalactic 



GMC samples (e.g. Solomon et al. 1987 : Bolatto et al. 2008h . However the MAGMA LMC clouds are offset towards narrower 



linewidths and lower CO luminosities compared to GMCs of a similar size in these samples. The scatter in the scaling relations 
corresponds to order of magnitude peak-to-peak variations in the CO-to-H2 conversion factor (as inferred from the ratio of 
the virial mass to the CO luminosity), the H2 mass surface density and the CO surface brightness of the MAGMA GMCs. 



3. The physical properties of star-forming GMCs are very similar to the properties of GMCs without signs of massive star 
formation. Sightlines through non-star-forming GMCs tend to have lower peak CO brightness, suggesting that the filling 
fraction and/or brightness temperature of the CO-emitting substructure is lower for clouds without star formation. 

4. We find a significant positive correlation between the peak CO brightness and CO surface brightness of the MAGMA 
clouds and the stellar mass surface density. We propose that these correlations are due to an increase in the CO brightness 
temperature and/or an increase in the abundance of CO relative to H2 in the stellar bar region. 

5. The velocity dispersion (ctv) of the MAGMA GMCs increases in regions with high Hi column density (A^(Hi)). Higher vol- 
ume densities and/or higher virial parameters for GMCs in regions with high 7V(Hi) could produce the observed correlation, 
although the MAGMA data does not provide unambiguous evidence for either of these alternatives. 



6. There is some evidence that the H2 mass surface density of the MAGMA LMC cl ouds increases with the interstellar kinetic 
pressure, Pext- Although the molecular cloud model proposed by ( Elmegreen 19891 ) predicts a relation between P^xt and the 
mass surface density of an atomic+molecular cloud complex, the MAGMA clouds do not fulfil the predictions of the model 
for reasonable values of the metallicity, radiation field and GMC envelope mass. 
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